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ABSTRACT 


This  paper  presents  the  results  of  an  exhaustive  search  to 

find  optimal  full  period  multipliers  for  the  multiplicative  con-  „  , 

f>,  *'*■  f'  ‘ 

■31 

gruential  random  number  generator  with  prime  modulus  2  -1. 

Here  a  multiplier  is  said  to  be  optimal  if  the  distance  between 
adjacent  parallel  hyperplanes  on  which  k-tuples  lie  does  not  ex¬ 
ceed  the  minimal  achievable  distance  by  more  than  25  percent  for 

k=2 . 6.  This  criterion  is  considerably  more  stringent  than 

prevailing  standards  of  acceptability  and  leads  to  a  total  of 
only  414  multipliers  among  the  more  than  534  million  candidate 
multipliers . 

Section  1  reviews  the  basic  properties  of  linear  congruential 

generators  and  Section  2  describes  worst  case  performance  mea 

<► 

sures.  These  include  the  maximal  distance  between  adjacent  par 
allel  hyperplanes,  the  minimal  number  of  parallel  hyperplanes, 
the  minimal  distance  between  k  tuples,  the  lattice  ratio  and  the 
discrepancy.  Section  3  presents  the  five  best  multipliers  and 


compares  their  performances  with  those  of  three  commonly  employed 
multipliers  for  all  measures  but  the  lattice  test.  Comparisons 
using  packing  measures  in  the  space  of  k-tuples  and  in  the  dual 


9  r  ’ 

3  l* 

V 


space  are  also  made.  Section  4  presents  the  results  of  applying 
a  battery  of  statistical  tests  to  the  best  five  to  detect  local 
departures  from  randomness.  None  were  found.  The  Appendix  con¬ 
tains  a  list  of  all  optimal  multipliers.  - 

Kev  words:  Congruential  generator,  discrepancy,  lattice  test, 


random  number  generation,  spectral  test. 
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INTRODUCTION 


This  paper  presents  the  results  of  an  exhaustive  search  to 

find  optimal  multipliers  A  for  the  multiplicative  congruential 

random  number  generator  Z^  =  A  Z^  ^  (mod  M)  with  prime  modulus 
31 

M  =  2  -  1.  Since  Marsaglia  (1968)  showed  that  k- tuples  from 

this  and  the  more  general  class  of  linear  congruential  generators 
lie  on  sets  of  parallel  hyperplanes  it  has  become  common  practice 
to  evaluate  multipliers  in  terms  of  their  induced  hyperplane 
structures.  This  study  continues  the  practice  and  regards  a 

multiplier  as  optimal  if  for  k  =  2 . 6  and  each  set  of  parallel 

hyperplanes  the  Euclidean  distance  between  adjacent  hyperplanes 
does  not  exceed  the  minimal  achievable  distance  by  more  than  25 
percent.  The  concept  of  using  this  distance  measure  to  evaluate 
multipliers  originated  in  the  spectral  test  of  Coveyou  and 
NacPherson  (1967)  and  has  been  used  notably  by  Knuth  (1981). 
However,  the  criterion  of  optimality  defined  here  is  considerably 
more  stringent  than  the  criteria  that  these  writers  proposed.  In 
fact,  among  the  more  than  534  million  full  period  multipliers  A 
examined  in  this  study,  our  search  identified  only  414  optimal 
multipliers. 

First  proposed  by  Lehmer  (1951),  the  multiplicative  congruen¬ 
tial  random  number  generator  has  come  to  be  the  most  commonly  em¬ 
ployed  mechanism  for  generating  random  numbers.  Jannson  (1966) 
collected  the  then  known  properties  of  these  generators.  Shortly 
thereafter  Marsaglia  (1968)  showed  that  all  such  generators  share 


a  common  theoretical  flaw  and  Coveyou  and  MacPherson  (1967), 
Beyer,  Roof  and  Williamson  <1971),  Marsaglia  (1972)  and  Smith 
(1971)  proposed  alternative  procedures  for  rating  the  seriousness 
of  this  flaw  for  individual  multipliers.  Later  Niederreiter 
(1976,  1977,  1978a, b)  proposed  a  rating  system  based  on  the  con¬ 
cept  of  discrepancy,  a  measure  of  error  used  in  numerical  inte¬ 
gration.  With  regard  to  empirical  evaluation,  Fishman  and  Moore 
(1982)  described  a  comprehensive  battery  of  statistical  tests  and 
illustrated  how  they  could  be  used  to  detect  local  departures 
from  randomness  in  samples  of  moderate  size  taken  from  these 
generators . 

Although  the  theoretical  rating  procedures  have  existed  for 

some  time,  with  the  exception  of  Hoaglin  (1976),  Ahrens  and 

Dieter  (1977)  and  Knuth  (1981),  little  use  has  been  made  of 

them.  The  present  study,  by  its  sheer  exhaustiveness,  removes 
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this  deficiency  for  generators  with  M  =  2  -  1.  Section  1 

reviews  the  basic  properties  of  linear  congruential  generators. 
Then  Section  2  describes  the  worst  case  performance  measures  that 
have  been  proposed  to  rate  generators  in  k  dimensions.  These 
include  the  maximal  distance  between  adjacent  parallel  hyper¬ 
planes,  the  minimal  number  of  parallel  hyperplanes,  the  minimal 
distance  between  k- tuples,  the  lattice  ratio  and  the  discre 
pancy.  These  concepts  are  described  in  this  study  principally  in 
terms  of  the  space  of  k- tuples  and,  where  appropriate,  in  terms 
of  the  dual  lattice  space.  However,  in  order  not  to  obfuscate 
central  concepts  the  exposition  relies  on  a  minimal  use  of  formal 
lattice  theory. 
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Section  3  presents  the  five  best  multipliers  and  compares 
their  performances  with  those  of  three  commonly  employed  multi¬ 
pliers  for  all  these  measures  but  the  lattice  test.  The  Appendix 
contains  a  list  of  all  optimal  multipliers.  Also,  lattice  pack¬ 
ing  measures  are  presented  and  again  show  the  dominance  of  the 
five  best  over  the  three  commonly  used  multipliers.  Packing 
measures  in  the  dual  space  are  also  computed.  This  last  concept 
is  identical  with  Knuth's  figure  of  merit  for  evaluating 
generators.  Our  results  indicate  that  with  regard  to  this 
criterion  the  five  best  perform  better  than  all  30  multipliers 
listed  in  Table  1  of  Knuth  (1981,  pp.  102-103).  Bounds  on 
discrepancy  are  also  computed  and  discussed. 

Section  4  presents  the  results  of  a  comprehensive  empirical 
analysis  of  the  local  sampling  properties  of  the  best  five,  using 
the  procedures  in  Fishman  and  Moore  (1982).  No  evidence  of 
departures  from  randomness  was  detected. 

1.  Linear  Conuruential  Generators 

A  linear  congruential  generator  produces  a  sequence  of  non¬ 
negative  integers 

{ZQ,  Z.  =  AZ.  +  C  (mod  M);  i=l,2,...}  (1) 

where  the  modulus  M,  and  multiplier  A  are  positive  integers  and 
the  seed  ZQ  and  constant  C  are  nonnegative  integers.  For  pur¬ 
poses  of  conducting  sampling  experiments  on  a  computer,  the  ele¬ 
ments  of  the  sequence  Z  are  normalized  to  produce  the  sequence 

U  *  (Oj  s  Z^/M;  i«l,2, . . . }  <2> 
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whose  elements  are  treated  as  if  they  were  sampled  independently 
from  the  uniform  distribution  on  the  interval  [0,1).  The  objec¬ 
tive  in  assigning  values  to  M,  A,  ZQ  and  C  is  to  make  the 
errors  incurred  in  this  treatment  of  U  tolerable  ones.  Here 
errors  are  principally  of  two  types,  one  being  the  approximation 
of  a  continuous  phenomenon  on  (0,1)  by  the  discrete  sequence  U 
and  the  other  being  the  distributional  distortions  in  U  induced 
by  the  use  of  the  deterministic  generator  (1).  In  addition,  com¬ 
putational  considerations  play  a  role  in  choosing  M,  A  and  C. 

One  property  of  the  generator  (1)  is  the  period 

T  =  min  {k>l:  Z  .  =  7  for  all  n  >  M} .  (3) 

n+k  n 

The  larger  M  is,  the  larger  T  can  potentially  be,  and  conse¬ 
quently  the  denser  the  points  of  U  are  in  [0,1).  The  more  dense 
these  points  are,  the  smaller  the  continuity  error  is. 

Table  1  lists  several  types  of  linear  congruential  generators 
that  are  or  have  been  in  common  use.  Here  A,  C,  Z^  in  the  table 

Insert  Table  1  about  here 

guarantee  maximal  period  for  the  corresponding  modulus  M.  Note 

that  types  1  and  2  give  full  periods  whereas  the  remaining 

generators  give  only  one  fourth  of  the  numbers  between  1  and 

2°.  Moreover,  types  4a  and  4b  do  not  produce  equidistributed 

B 

sequences.  Also,  the  use  of  M=2  enables  one  to  replace  divi¬ 
sion  and  multiplication  by  less  time  consuming  shift  and  add 

0 

operations.  Although  M=2  -1  does  not  allow  this  substitution 


/ 


I*. 


'Ll. 


**+**«+ 
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i: 


directly,  a  procedure  due  to  Payne,  Rabung  and  Bogyo  (1969) 

enables  one  to  retain  part  of  this  improved  efficiency.  Note 

M- 1  - 

that  A  is  a  primitive  root  of  M  if  A  =1  (mod  M)  and 
AQ  ^  1  (mod  M>  for  0<Q<M-1. 

Today  only  linear  congruential  generators  of  types  2  and 

3  are  commonly  used.  On  18M  computers  with  a  word  size  of 

32  bits  and  C  =  0,  the  generator  called  SUPKR-DUPER 

(Marsaglia  1972)  uses  M  =  232,  A  =  69069  and  ZQ  =  odd 

30 

integer  to  give  a  period  T  =  2  For  generators  of  type 
2  with  prime  number  modulus  M  =  231  -  1,  APL  (Katzan  1971) 
uses  A  =  16807,  the  S1MSCR1PT  II  programming  language 
(Kiviat,  Villanueva  and  Markowitz  1969)  uses  A  =  630360016, 
SAS  (1982)  uses  A  -  397204094  and  the  IHSL  Library  (1980) 

gives  the  user  the  choice  of  A  =  16807  or  A  =  397204094. 

31 

The  resulting  period  is  T  =  2  -2. 


2.  Theoretical  Measures  Of  Performance 


In  practice,  it  is  relatively  common  to  use  the  pseudorandom 
numbers  produced  by  (1)  in  groups  or  k-tuples .  Consider  the 
sequence  of  points 


W.  =  {V.  .  =  (Z.  . . Z.  .); 

k  i,k  l+l  l+k 

and  the  normalized  sequence 


i  =  1,2,...} 


(4) 


Vk=  lVi,k=  (Zi+l/M . Zi+k/M)i  1  <5> 

Ideally  one  wants  the  sequence  of  points  to  be  equidistri- 

buted  in  the  k-dimensional  unit  hypercube  for  k  =  2,3,....  How¬ 
ever,  the  form  of  the  generator  (1)  limits  the  extent  to  which 
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one  can  achieve  this  ideal.  For  example,  observe  that  an  ideal 
generator  of  the  integers  1  =  {1.....K  1}  produces  (M-l) 
equidistributed  points  in  k- dimensional  space  whereas  a  generator 
of  type  2  produces  only  M-l  points  in  the  this  space. 

Although  this  constancy  of  the  number  of  points  is  itself 
sobering,  it  is  one  of  two  important  issues.  To  illustrate  the 
second  issue,  Fig.  1  shows  a  plot  of  2-tuples  for  the  genera¬ 
tors  Z.  5  7  Z.  ,  (mod  61)  and  Z.  =  31  2.  ,  (mod  61)  where  61 
l  l-l  l  l-l 

is  a  prime  number  and  7  and  31  are  primitive  roots  of  61.  Al¬ 
though  no  one  would  seriously  use  either  of  these  generators  to 

Insert  Fig.  1  about  here. 

produce  random  2-tuples,  a  comparison  of  Figs,  la  and  lb  arouses 
a  concern  that  holds  for  more  realistic  generators  as  well. 

Notice  that  the  distribution  of  points  in  Fig.  lb  is  considerably 
less  uniform  than  the  distribution  in  Fig.  la.  Since  such  dif¬ 
ferences  in  two  and  higher  dimensions  are  attributable  entirely 
to  the  cl  of  multiplier  and  since  there  are  an  enormous 

number  of  candidate  multipliers,  a  deep  analysis  of  k-tuples 
generated  by  (1)  across  all  those  multipliers  is  needed  to  assess 
the  extent  to  which  the  resulting  sequences  depart  from  the 
ideal  of  equidistributedness . 

Several  theoretical  procedures  have  been  proposed  to  make 
this  assessment.  They  include: 

1.  maximal  distance  between  adjacent  parallel  hyperplanes 


(spectral  test) 
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2.  minimal  number  of  parallel  hyperplanes 

3.  minimal  distance  between  points 

4.  ratio  of  lengths  of  longest  and  shortest  minimal  basis 
vectors  (lattice  test) 

5.  discrepancy. 

Although  diverse  in  what  they  measure,  the  procedures  share  a 
common  unifying  concept.  All  follow  from  recognizing  that,  with 
the  exception  of  generators  of  types  4a  and  4b,  the  k-tuples 
can  be  regarded  as  points  in  a  regular  lattice.  Moreover, 
generators  of  types  4a  and  4b  lead  to  k-tuples  on  two  intermeshed 
regular  lattices.  Ahrens  and  Dieter  (1977),  Beyer,  Roof  and 
Williamson  (1971)  and  Coveyou  (1970)  provide  detailed  descrip¬ 
tions  of  this  relationship  with  lattice  theory.  To  keep  the 
focus  of  attention  on  the  assessment  of  interest,  the  present 
paper  presents  only  the  features  of  lattice  theory  that  are 
essential  for  describing  these  procedures.  Also,  unless 
otherwise  noted  our  description  applies  for  generators  of  type  2. 
Comparable  analyses  can  be  performed  for  each  other  type  of 
generator. 

2.1  Maximal  Distance  Between  Adjacent  Parallel  Hyperplanes 
Observe  that  Z.  can  be  written  in  the  form 

l 

Z.  =  Z^A1  (mod  M) 

l  0 


i-1 


M  I 

m=0 


i>l 


(6) 


i 


T-* 


where  =  |^AZ  ^/H  j=l  ,2,...  Now  f or  k  >  1 ,  q  =  (qQ, 


•v^ 


and  y  consider  the  k-  dimens ional  hyperplane 


HJc(q.  y)  -  {<*0 . xic_i):  E  “Jj*!  =  y) 

~  j=0  J  J 


and  in  particular  the  family  of  parallel  hyperplanes 


Hfc(q)  =  { < q ,  y) :  y  -  0  (mod  1)}. 


Observe  that  the  elements  of  V  in  (5)  lie  on  hyperplanes  in  H  (q) 

JC 


in  (7)  if 


<i>  qQ . qk  _1  integer 


(ii)  q(A)  =  Z  q.AJ  i  0  (mod  M)  . 


These  restrictions  are  sufficient  since  for  any  V.  .  in  V  the 

l  ,k  k 


quantity 


,  k-1  Z  A  ... 

y.  =  —  £  q  7  =  -9—  q(A)  +  ki 

yi  M  qj  i+d  M 


where 


k-1  i+j-1 


c.  =  -Z  q.  Z  K.  .  A‘  . 
1  j-1  3  m=0  1+Jm 


Restriction  (i)  insures  that  k^  is  an  integer  and  restriction 


(ii)  insures  that  y.-k,  is  an  integer.  These  restrictions 
i  l 


hold  throughout  the  remainder  of  this  paper. 


For  the  ensuing  analysis  it  is  convenient  to  extend  V. 


modulo  one  to  the  set 


Vfc  =  {V  =  (vQ, . . .  ,vk  l>  integer}  U  {V  =  (vQ . vjc-i)  -  vt  k<"*od  l>i 


i=l . T} 


1 


I 


Since  (i)  and  (ii)  hold,  the  points  in  v£  also  lie  on  hyperplanes 
in  H^Cq).  The  set  of  all  hypcrplanes  containing  at  least  one  point 
a 

of  is  then 

*  fc~l  ^  a  a 

H.(q)  =  {H  (q):  y  =  I  q.v.,  v  cV  } .  (11) 

~  ~  j*o  J  J 

Moreover,  one  can  index  these  hypcrplanes  by  the  set  of  integers 

k-i 

a  ,  a  a  a  a 

\<q)  -  {y  =  I  q.v  :  V  cVj  .  (12) 

j=0  J  J 

We  now  use  these  representations  to  show  that  for  specified 
k  >  1  and  q  the  k-tuples  in  lie  on  a  sot  of  parallel  hyperplanes 

for  which  the  Euclidean  distance  between  adjacent  hyperplanes  is 

a  a 

fixed.  The  set  of  hyperplancs  is  Hfc(q)  and  for  y  and  z  in  Yfc(q)  the 

k-1 

2  1/2 

Euclidean  distance  between  Hfc(q,  y)  and  H^(q,  is  |y-z|/(  I  q.) 

j=*0  3 

a 

To  prove  the  result,  it  suffices  to  show  that  the  Yk(q)  is  composed  of 
all  integer  multiples  of  some  fixed  constant  I^(q),  for  then  the 

a 

Euclidean  distance  between  adjacent  byperplanes  in  H^(q)  is 


dk<q;  A,  M) 


V3> 


k-1 

(  r 

j-0 


2.1/2 

V 


(13) 


a 

By  way  of  proof,  note  th«t  if  V  and  V'  are  two  elements  of  then 

a  a 

V* '  =  V*  -  V  is  also  •'  >  1  therefore  for  y  and  z  in  Y^(q)  one  has 


I 
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*  *  * 
y- z  m  Y^(q).  Also,  for  any  integer  j  and  point  V  in  the  point 

.  *  *  *  * 

V  =  jV  is  also  in  so  that  z=jy,  for  ycY^Cq),  is  also  in  Y^(q). 

* 

Therefore,  it  follows  that  all  elements  of  Y^(q)  are  multiples  of 

I^q)  «  min{  I y* | >0;  y*cY*(q>), 

thus  establishing  (13).  Without  loss  generality  we  take 
(iii)  1  (q)  =  1. 

Since  many  different  vectors  q  satisfy  (i),  (ii),  and  (iii) 

for  a  given  multiplier  A  and  induce  families  of  parallel  hyper- 
planes,  additional  criteria  are  needed  to  enable  one  to  charac¬ 
terize  the  extent  of  equidistributedness  of  the  k- tuples  in 
(5)  in  the  k-dimensional  unit  hypercube  for  each  possible  multi¬ 
plier.  One  such  criterion  is  the  maximal  distance  between  adja 
cent  parallel  hyperplanes  which  is  a  worst  case  measure  for  a 
particular  multiplier  A.  It  is 

k-1  ,  -1/2 

d*  (A,  M)  =  max  1<X  qT>  1  (14) 

k  qo . Vl  j  -0  J 

subject  to  restrictions  (i),  (ii)  and  (iii).  In  particular,  note 
that  the  constraint  (iii)  eliminates  the  numerator  of  (13)  from 
the  maximization  (14). 

When  using  (14)  to  compare  k-tuple  performance  for  several 
multipliers  for  a  type  of  generator,  one  prefers  the  multiplier 


that  gives  the  minimal  maximal  distance  since  this  implies 


smaller  empty  regions  in  the  k- dimens ional  unit  hypercube  for 
this  multiplier  than  for  the  other  multipliers.  However,  there 
is  a  limit  to  how  small  this  maximal  distance  can  be;  in 
particular,  it  is  known  that  (Cassels  1959,  p.  332) 


Ml/1Cd*(A,M)  >  =  •!  2  k=4  (15 


(3/4)1M 

k=2 

2-l/6 

k=3 

2-1M 

k=4 

2-3/10 

k=5 

(3/64)1/12 

k=6 

To  illustrate  the  significance  of  these  bounds,  note  that 
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with  the  modulus  M  =  2  -  1  one  has 


*  31 

d.  (A,  2Ji 


1)  > 


2008 

X 
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i 

o 

k=2 

6905 

X 

ro 

l 

o 

k=3 

3906 

X 

CM 

o 

r-t 

k=4 

1105 

X 

10  1 

k=5 

2157 

X 

10-1 

k=6 , 

indicating  the  relative  coarseness  of  the  grid  of  points  in  as 
few  as  four  dimensions. 

Using  multivariable  Fourier  analysis,  Coveyou  and  McPherson 
(1967)  advocated  using  the  minimized  "wave  number" 


k-1  ,  1/2 
<  l  q:> 
j=0  3 

(s.t.  q  , . . . ,q^  ^  integer  and  q(A)  =  0  (mod  M) )  to  determine 
the  relative  desirabilities  of  alternative  multipliers,  hence  the 
name  spectral  test.  Shortly  thereafter,  it  became  apparent 
(Coveyou  1970;  Beyer,  Roof  and  Williamson  1971)  that  one  could 
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perform  equivalent  studies  using  (14)  by  viewing  the  k-tuples  as 
being  arranged  on  parallel  hyperplanes  and  exploiting  the  mathe 
matical  properties  of  the  so  induced  lattice  structure .  In  fact, 
it  turns  out  that  the  physical  interpretation  of  results  can  be 
more  easily  understood  in  the  space  of  whereas  the  computa 
tional  procedures  are  more  easily  understood  by  working  in 
the  dual  space  of  q.  We  return  to  this  issue  in  Section  3. 


2.2  Minimal  Number  of  Parallel  Hyperplanes 

A  second  measure  of  equidistr ibutedness ,  suggested  by 
Marsaglia  (1968) ,  is  the  number  of  parallel  hyperplanes 

Nk^q0’ ‘ ‘ ,qk  l*  A’  on  w*l*c*1  the  ^-tuples  lie.  If 
this  number  is  small  for  a  particular  multiplier  A,  then  this  is 
an  indication  that  there  exist  large  regions  in  the  k-dimensional 
unit  hypercube  that  contain  no  k-tuples. 

Observe  that  with  restriction  (iii)  gives  the  upper  bound 

k-1 

Nfc(q0 . qk_r  A,  M)  <^I  |qj|.  (16) 


Using  the  development  in  Dieter  (1975),  one  also  observes  that 
k-1  _  k-1 

-  I  (9.)  <  <  l  (qj  +  U1.....T 

j=0  J  1  j-0  3 

where 

10  if  x  >  0 

-x  if  x  <  0 

and 

I  x  if  x  >  0 


+ 

x  « 


0  if  x  <  0. 


f 
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Because  of  the  restrictions  (i)  through  (iii)  the  number  of  dis¬ 
tinct  Yj  (mod  1)  is  precisely  the  number  of  parallel  hyperplanes 
Nk<<J0 . qk_i;  A*  M><  then 

k- 1  k- 1  k- 1 

Vqo . q.  *.  A,  M)  =  l  (q.)~  +  l  <q.)  +  =  £  |q.l  -1.  (17) 

*  0  j=0  3  j=0  3  j  -0  3 

The  -1  in  (17)  follows  from  the  observation  that  Y.  takes  on 

l 

all  integer  values  (qQ . qk-l*  A'  • 

As  before,  there  exist  many  vectors  q  that  satisfy  restric¬ 
tions  (i),  (ii),  and  (iii).  A  worst  case  measure  here  is 

N*(A,  M)  =  min  ■  ■  ■  .qfc ^  A,  M)  (18) 

subject  to  the  four  restrictions.  When  using  this  criterion  to 

choose  among  several  multipliers,  one  prefers  the  one  that  gives 

the  maximal  minimal  '  '  ‘  ,qk-l*  '  As  *n  the  case 

of  distance  between  hyperplanes,  an  upper  bound  exists  on 
* 

N^A,  M) ,  namely  (Marsaglia  1968) 


*  1/k 

Nk<A,  M)  <  (k!M>  k=l ,2 . 

31 

In  particular,  for  M  =  2  -  1  the  bounds  are 

1  65536  k -2 

2344  k=3 

476  k=4 

191  k=5 

107  k=6 . 


Again,  these  bounds  are  limiting,  and  encourage  one  to  search  for 
multipliers  that  can  come  close  to  the  bounds. 

Knuth  (1981,  p.  92)  points  out  that  the  ordering  of  several 


multiplies  A^,...A  according  to  the  maximal  distance  measure 


* 

d^(A,M)  may  differ  from  the  ordering  established  by  the 

* 

minimal  number  of  parallel  hyperplanes  measure  N^(A,M).  In 

* 

particular,  he  notes  that  N^(A,M)  "is  biased  by  how  nearly 

the  slope  of  the  lines  or  hyperplanes  matches  the  coordinate  axes 

* 

of  the  cube."  That  is,  Nfe(A,M)  may  be  relatively  large  when 

ft 

dfc(A,M)  is  also  relatively  large.  Since  in  this  case  one 

inclines  to  discount  the  multiplier  because  of  sparseness  indi 
* 

cated  by  d^lA.M) ,  there  is  some  justification  for  valuing 
*  * 

d^(A,M)  more  highly  than  Nfc(A,M)  as  a  measure  of  per¬ 
formance.  Section  3  takes  this  into  consideration  when  searching 

for  optimal  multipliers. 

* 

Although  d^iA.M)  provides  a  more  definite  evaluation  of 
a  multiplier  than  Nfc(A,M)  does,  the  latter  quantity  has  at 
least  one  readily  appealing  attribute  that  justifies  its  con¬ 
sideration.  We  illustrate  this  feature  for  the  type  3  generator 
with  A=65539  and  M=231.  This  generator  is  known  as  RANDU  and 
was  a  standard  feature  of  the  IBM  Scientific  Subroutine  Library 
on  360/370  series  computers  for  many  years.  Observe  that  65539  = 
216+3  so  that 

Z.  ,  =  (216  +  3)Z.  <mod  231) 

1+1  1 

16  31 

Z.  -=(6x2  +  9>Z.  (mod  2  ) 

i+2  i 

31 

Z.  „  =  6Z.  ,  -  9Z.  (mod  2  ) 

i+2  i+l  i 

31 

Z.  „  -6Z.  ,  +  9Z.  =  0  (mod  2'  > . 

i+2  i+l  i 

Moreover, 

U.  „  -  6U.  ,  +  9U.  r  0  (mod  1) 

i+2  i+l  i 

indicating  that  N*(65539,  231)  <  16,  a  devastating  indictment 


of  RANDU  in  three  dimensions.  Therefore,  the  valuable  feature  of 
* 

is  that  it  can  on  occasion  identify  a  poor  multiplier  with 
relatively  little  computational  effort. 

2.3  Distance  Between  Points 

Smith  (1971)  has  suggested  an  alternative  measure  of  equidis 
tributedness  based  on  the  minimal  distance  between  k- tuples 

k-1 

*  l  2  1/2 

c  (A,  M)  =  min  [  i  l  (2.  .  -  Z  .)  ]  .  (19) 

k  n  .  _  M  .  l+j  m+j 

l<i,m<T  j  =0 

i  m 

Since  the  total  number  of  points  is  fixed  at  T,  the  smaller 
* 

c|t(A,  H)  is  for  a  given  A,  the  larger  the  volume  of  space  is 
in  the  k-dimensional  hypercube  that  is  devoid  of  points.  There¬ 
fore,  when  comparing  several  multipliers  in  k  dimensions  one  pre- 

* 

fers  the  one  that  gives  the  maximal  c^tA,  M) . 

* 

Whereas  d^(A,M)  measures  distance  between  adjacent 
* 

parallel  hyperplanes,  ck(A,  M)  measures  distance  between 

nearest  points.  An  alternative,  but  equivalent  interpretation  is 
* 

to  view  l/c^(A,  M)  as  the  maximal  distance  between  adjacent 
parallel  hyperplanes  in  the  dual  space  of  <j.  This  observation 

A 

enables  one  to  establish  the  upper  bounds  for  ck(A,  M) 
k=2,3 . . .  (Cassels  1959,  p.  332): 

*  1/k 

cfc  (A,  M)  <  l/yfc  M  (20) 

where  is  defined  in  (15).  This  duality  relationship  also 

A 

facilitates  the  computation  of  ck<a,  M)  using  the  algorithm 


i 

, 


in  Dieter  (1975) . 
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2.4  Discrepancy 

The  concept  of  discrepancy  originated  in  the  study  of  how 
well  equidi stributed  sequences  perform  in  sampling  procedures 
designed  to  approximate  the  volumes  of  regions  in  the  k-dimen 
sional  unit  hypercube  and  in  numerical  integration.  Having 
recognized  the  relationship  between  this  problem  and  that  of 
measuring  the  performance  of  a  random  number  generator, 
Niederreiter  (197/)  adapted  the  discrepancy  measure  to  this 
latter  problem  and  gave  bounds  for  it. 

Consider  the  sequence  of  k-tuples  {W.  ;  i —1 ..... T }  de- 

l ,  x 

fined  in  (4).  For  N-1.....T  discrepancy  in  fc  dimensions  for  a 
multiplier  A  and  a  modulus  M  is  defined  as 


D^k>(A,M) 


max 

R 


number  of  tf  . wu  R 

N 


volume  of  R  1 
Mk 


where  R  ranges  over  all  sets  of  points  of  the  form  R  = 

{(w_ , . . .  ,w.)  | a  <  w  <  8  , . . .  ,a  <  w.  <  8.  } .  Here  a.  and 
1  k  l  l  1  k  k  k  j 

8.  are  integers  in  the  range  0  <  a.  <  B,  <  M  for  1  <  j  <  k  so 
J  1  j 

that  the  volume  of  R  is 


k 

,n. (8.  -  a.). 

j=l  J  3 

Niederreiter  (1977,  1978a)  gave  upper  and  lower  bounds  for 
(k) 

D  (A,H)  for  generators  of  types  1,  2  and  3  for  arbitrary 

N 

N  <  T.  In  particular,  the  upper  bound  for  generators  of  type 
is 


\ 


'  w  < 


I 


(21) 


D'  ( A,M)  <  k  +  min(N.  (M-W)  )  I  _ 

M  H  ,  .  .  r(q,M) 

q  (mod  m)  _ 


1/2 

+  max  (0.  N-  (M-N)1'*)  z* 


q  (mod  M) 
q(A)  Ir  0  (mod  M) 


r(q ,  M) 


where  the  asterisk  denotes  exclusion  of  q^=  . . . =q^  ^  -0, 


k-1 


r(g,  M)  =  I  r(q  . ,  M) , 
j  =0  J 


r(q,M)  *  1 


if  q  r  0  (mod  M) 


-  M  sinfllq/MlI  if  q/3  0  (mod  M) 

and 

II til  »  min  (t,  1-t) . 

Note  that  this  bound  holds  for  any  local  sample  of  N  successive  k  tuples  from 
the  generator  as  well  as  for  a  global  evaluation  of  performance  when  N=T. 

At  present  there  exists  no  algorithm,  other  than  total  enu¬ 
meration,  for  computing  this  upper  bound  and  this  situation  is 
likely  to  remain  so.  However,  the  form  of  the  bound  enables  one 
to  establish  a  valuable  relationship  between  the  spectral  test 
and  discrepancy.  Note  that  sinllq/Mll>0  in  the  bound  and  that 
the  number  of  such  terms  is  a  function  of  q(A) .  Recall  that  the 


k-1 

(  X 


2 


,1/2 


quantity 


r 


-18- 


is  minimized,  subject  to  q(A)  =  0  (mod  M) ,  to  find  the  maximal 
distance  between  adjacent  parallel  hyperplanes.  If  this  mini¬ 
mized  quantity  turns  out  to  be  small  for  a  multiplier  A,  then  the 
congruence  occurs  frequently  for  0  <  |q^  ^ |  <  M  j=l,...,k. 

This  clearly  adds  positive  terms  to  the  second  summation  and  there 
fore  the  upper  bound  is  large.  If  for  an  alternative  multiplier  A' 
the  minimized  quantity  turns  out  large,  the  congruence  holds  less 
frequently  and  the  upper  bound  is  smaller  than  in  the  previous  case. 
Thus  the  results  for  the  spectral  test  convey  useful  information 
about  the  bounds  on  discrepancy. 

For  generators  of  types  1  and  2,  Niederreiter  (1976,  1978a)  also 
gave  the  lower  bound 

D<k><A.K>  ,J  I'kV^A.M  for  2  <  k  <6 

(  */2<2i r  +  1) V  '(A.M)  for  k  >  7 

and  the  upper  bound 


dk)..  k  .  ...  /M-T  2  _  „  7.k 

Dt  (A.M)  <  -  +  min  (1,  )(-  log  M  +  -) 


1-k  k  k-1  (k) 

+  (log  2)  ((2  log  M)  +  4  <2  log  M)  )/2p  (A.M) 


(23) 


_k  ,_k-2  , ,  .J+k-2..  (k)..  ... 
+2  (2  -1)  (  >/p  (A.M) 


where 


p(lc)(A,  M)  = 


min 

q(mod  M) 
q~  *  (0 . 0) 

q(A)  r  0  (mod  M) 


k-1 

(  n  max(l ,  lq  .  I ) )  . 
j=0  J 


(24) 


and 


(log  M) /log  2 
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Comparable  results  exist  for  generators  of  type  3. 

With  the  exception  of  k=2  no  algorithm  exists  for  computing 
(k) 

p  (A,H) .  Ahrens  and  Dieter  (1977,  Theorem  5.17)  gave  the 
stronger  lower  bound 


D*k)(A,M)  >  l/£  min  (Xm  n  |q.  |)j 

q  *  (0 . 0)  1=0 

q(A)  =  0  mod  M 

where  m  denotes  the  number  of  non  zero  qj, 

m  ,  .  . 

m  if  m  =  2  or  3 

mm/(m-l)m  H  if  m  >  4 
ID  ~ 

and 

Lm/2j  +1 

Hm  *  ('1)J  (j>(Lm/2J  +l-j> l/(m-l)!l"  . 

For  k=2  Borosh  and  Niederreiter  (1983)  showed  that 


(25) 


(26) 


p(2*( A,M)  =  min  (|q  |*|tM-  q  A|) 

o  <  Iq-L !  <  M/2 

for  some  t  satisfying  qQ  =  tM  -  q^A.  This  result  makes  the 
bounds  in  (22)  and  (23)  operative  for  k=2. 

Niederreiter  (1977,  1978b)  provided  additional  bounds 
for  k=2.  For  type  2  generators 

(2)  P 

'(A,M)  <  (2  +  l  a^/T 
i=l 

and 


(27) 


(28) 


Dj  (A,M)  <  [2  +  C(K>  log  T]/T 


(29) 


-20- 


where  are  the  partial  quotients  in  the  conti 


nued 


fraction  expansion  of  A/M,  K  -  max  (a . ,a  )  and  C(K)  = 

1  P 

2/log  2  for  1  <  K  <  3  and  C(K)  =  (K+l)/log(K+l)  for  K  >  4.  Ex¬ 
pressions  (28)  and  (29)  also  hold  for  type  3  generators  with  2/T 

replaced  by  1/T  and  with  a^-.-.a  being  the  partial  quo~ 

B-2 

tients  of  A/2  Earlier,  Dieter  (1971)  derived  closely 

related  results  based  on  continued  fractions  to  nearest  integers 
rather  than  regular  continued  fractions. 

Borosh  and  Niederreiter  (1983,  Table  2)  have  carried  out  a 
systematic  search  for  multipliers  of  type  3  and  type  4.  In  par¬ 
ticular,  they  gave  maximal  period  multipliers  with  K  <  3  for 
B  -  6, 7,..., 35  for  each  type. 


2.5  Lattice  Test 

Beyer,  Roof  and  Williamson  (1971)  and  Marsaglia  (1972)  pro¬ 
posed  an  alternative  figure  of  merit,  for  evaluating  alternative 
multipliers,  based  on  the  concept  of  squareness .  We  use  Fig.  1 
to  illustrate  this  concept.  Clearly  one  can  construct  a  vast 
number  of  parallelograms  of  varying  areas  that  include  no 
interior  points.  The  presumption  of  the  lattice  test  is  that  one 
prefers  multipliers  that  produce  parallelograms  of  minimal  area 
whose  sides  are  close,  if  not  equal,  in  length;  hence,  the  notion 
of  squareness,  where  angles  are  neglected. 

Now  the  minimal  volume  of  a  k-dimensional  parallelepiped 

k-1 

generated  by  k-tuples  from  (1)  subject  to  (ii)  is  M  .In 
evaluating  a  particular  multiplier,  the  objective  of  the  lattice 
test  is  to  find  the  basis  vectors  that  come  closest 


/ 
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in  k  dimensions  to  achieving  this  squareness  for  parallelepipeds 
k-1 

of  volume  M  To  measure  the  extent  of  the  departure  from 

equidistributedness  in  k  dimensions  Beyer,  et  al.  and  Marsaglia 

recommended  the  quantity 

max  |a. | 

1  <  i  <  k  1 

R  (A)  = - ;  (30) 

min  |a. | 

1  <  i  <  k  1 

that  is,  the  ratio  of  the  lengths  of  the  longest  and  the  shortest 
basis  vectors.  Clearly  R^(A)  >  1  and  presumably  one  prefers 
multipliers  for  which  R^(A)  is  close  to  unity. 

It  is  worthwhile  noting  that  the  basis  vectors  a  , . . . ,o.  play 
an  implicit  role  in  the  previously  mentioned  tests  as  well.  For  example, 
one  can  show  that  for  k=2  the  maximal  distance  between  parallel  hyper- 
planes  is 

1/2 

d*  <A,H)  =  (|a2|2  -  |ox  •  o2|2/|0l|2) 


where  we  take  to  be  the  longer  vector. 

Although  the  figure  of  merit  in  <30)  has  intuitive  appeal, 
there  is  no  universal  agreement  about  its  usefulness  in  identify¬ 
ing  good  multipliers.  Marsaglia  (1972,  p.  275)  suggested  a  gene- 

32 

rator  of  type  3  called  super-duper  with  M  =  ?  and  A  -  69069. 

It  has  R2<A)  =  1.06,  R3(A>  =  1.29,  R^(A)  =  1.30  and  RgU)  =  1.25; 
an  appealing  generator  as  evaluated  by  the  lattice  test.  For  this  generator 
Niederreiter  (1978,  pp.  1027-1028)  showed  that  p*2*(A,M)  <  69069 

so  that  (22)  gives  D^2)(A,M)  >  1/(4  x  69069)  =  .3620  x  10"5  . 


But 
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Borosh  and  Niederreiter  gave  a  multiplier  A  =  3039177861  for  M  =  2^2 

p 

(2) 

with  p  (A,M)  -  .2517M  and  ][  a^  -  51  for  which  (28)  based  on 

M  =  23?  gives  D^,2\a,M)  <  (1  4  51)/230  -  .4843  x  10  7  .  This  result 

illustrates  that  although  super  duper  has  the  appealing  figure  of 

merit  R^(69069)  1.06,  there  exist  multipliers  with  A  r  5 

32 

(mod  2  )  that  dominate  it  by  a  substantial  margin  in  k-2 

dimensions  with  regard  to  discrepancy. 

3 .  Analys i s 

This  section  presents  results  of  an  investigation  based  on 
* 

the  evaluation  of  {dfc  (A,M);  k=2 . 6}  for  all  multipliers  A 

31 

that  are  primitive  roots  of  M=2  1,  using  an  algorithm  of 

Dieter  (1975),  as  described  in  Knuth  (1981,  algorithm  S).  Hardy 
and  Wright  (1960)  show  that  the  number  of  primitive  roots  for  M 
prime  is  <J>(M-1)  where 

<J)(M-1)  =  number  of  integers  not  exceeding  and 
relatively  prime  to  M-l. 

This  quantity  is  called  the  Euler  totient  function.  Since 

4>(M- 1 ) / (M-l )  i  .249  for  M=231-l  (Ahrens  and  Dieter  1977,  p.  7,6), 

31 

one  has  4»(2  -2)  =  534723428,  a  not  inconsequential  number. 

To  find  the  primitive  roots,  one  notes  that  if  B  is  the 
smallest  primitive  root  of  the  prime  modulus  M,  then  every 
primitive  root  has  the  form 

A  r  B1  (mod  M) 


where  I  is  an  integer  whose  largest  common  factor  with  M- 1  is 

unity.  Since  one  also  can  show  that  for  every  such  I  there  exists 

1  MIX 

a  pair  of  multipliers  B  (mod  M>  and  B  (mod  M)  with  identi 

cal  lattice  structures,  it  suffices  to  investigate  only  half  of  all 

the  primitive  roots.  In  the  present  case  7  is  the  smallest  primi- 
31 

tive  root  of  2  -1  so  that  only  267361714  multipliers  require 

examination.  Note  that  the  multiplier  with  exponent  M- 1  1  produces 
the  same  sequence  as  the  multiplier  with  exponent  1  does,  but  in 
reverse  order. 

Clearly  one  needs  to  adopt  a  screening  procedure  to  identify 
and  collect  those  multipliers  that  "perform  well".  For  present 
purposes,  the  multipliers  of  most  interest  are  those  that  "perform 

well"  in  k.^2 . 6  dimensions  relative  to  the  constraints  imposed 

on  all  lattices  in  these  dimensions.  Consider  the  ratios 
Slilc(A,M)  =  Y|c/d£<A,M)Ml/k  k=2,...,6. 

As  seen  from  (15),  0  <  jc( A,M)  <  1.  Now  the  closer 
Si.i(A.M) . Si,6(A.M>  are  to  unity  the  better  the  perfor¬ 

mance  is  of  this  multiplier  with  regard  to  the  achievable  bounds 
in  2,..., 6  dimensions.  Therefore,  one  way  to  perform  the  screen¬ 
ing  is  to  identify  all  multipliers  for  which 

min  y(A,M)  >  S  0  <  S  <  1 

2<k<6 

where  S  is  specified. 

Initially  we  chose  S=.75.  Since  preliminary  computations 
indicated  that  there  were  an  unmanageable  number  of  multipliers 
that  satisfied  this  criterion,  we  changed  S  to  .80.  This  result¬ 
ed  in  a  total  of  207  optimal  multipliers,  as  listed  in  the 
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Appendix.  Recall  that  there  are  actually  twice  this  number  of 
optimal  multipliers.  The  abrupt  reduction  in  the  number  of  opti¬ 
mal  multipliers  when  shifting  from  S=.75  to  S=.8  is  itself  not¬ 
able.  Also  note  that  any  multiplier  for  which  fc(A,M)  >  .8 
for  k=2,...,6  guarantees  that  for  each  k  the  distance  between 
adjacent  hyperplanes  does  not  exceed  the  minimal  achievable 
distance  by  more  than  25  percent. 

Kor  each  selected  multiplier  and  k=2,...,6  we  also  computed 
the  ratios 

S2  fc(A,M)  -■  Nk(A,M)/(k!M>lc 

and 

S3  k(A,M)  =  ck(A,M)/YkM1/1C, 
again  using  Dieter’s  algorithm. 

Table  2  presents  these  ratios  for  the  multipliers  with  the 
five  largest  min  S  (A,M) .  It  also  presents  results  for 

1  v  K 

A=16807  which  is  in  APL  and  IMSL,  for  A=397204094  which  is  in 
1MSL  and  SAS,  for  A -630360016  which  is  the  SIMSCK1PT  11 
multiplier,  and  for  A=7.  This  last  multiplier  illustrates  the 
contrasts  that  are  possible  in  performance. 

Table  2  allows  one  to  make  several  notable  observations: 

Insert  Table  7  about  here. 

(a)  The  first  five  multipliers  perform  considerably  better 
than  the  remaining  multipliers  in  the  table  with  regard 
to  the  screening  measures  {S  ( A,M) }  and  with 

1  i  1C 

regard  to  {S?  k< A,M) }  and  {S3  ^(A.M) } . 


f 
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(b)  For  each  of  these  five  multipliers  S  (A,M) . 

*■  %  c 

S.  ,(A,M)  are  remarkably  close. 

1,0 

<c)  The  measures  S„  „(A,M) , . . . ,S,  , (A,M)  are  also  re 
3,2  3,o 

markably  close  and  behave  essentially  as 
S,  0(A,M) . . . ,S  (A,M) .  As  expected,  S  _(A,M)  = 

1)2  1  |  D  1  |  4 

s3  2<a,m>. 

(d)  S  (A,M),...,S  (A,M)  show  considerably  more 

2,7  2,6 

variation;  no  doubt  a  reflection  of  the  suboptimality  of 

these  multipliers  with  regard  to  this  criterion. 

We  now  turn  to  another  method  of  evaluating  performance  which 

derives  from  the  concept  of  packing  a  lattice  with  spheres  (see 

* 

Cassels  1959).  Recall  that  cfc(A,M)  is  the  distance  between 

nearest  points  in  the  unit  hypercube  of  k  tuples .  Then  the 

volume  of  a  sphere  with  this  diameter  is 

k/7  *  k 

ir  (c  (A,M>/2) 

Lk<A,M)  =  r<k/2  7T) 

where  T(»)  denotee  the  Gamma  function.  Suppose  one  packs  the 
lattice  with  such  spheres  centered  on  each  of  the  M  1  points 
in  (5)  and  at  the  origin.  Note  that  these  spheres  merely  touch 
and  that  6ince  there  are  only  M  k- tuples,  the  proportion  of  the 
unit  hypercube  packed  with  these  spheres  is  ML^tA.M). 

Let 

u^tA.M)  =  ?kMLk<A,M) . 

Using  the  lattice  packing  constants  in  (15)  and  (20)  one  has 
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)3 . 63  k=2 

5.92  k-3 

9.87  k^4 

14.89  k=5 

23.87  k=6  . 


Table  3  lists  u)^(AtM)  for  the  five  best  and  the  three 
other  commonly  employed  multipliers.  The  benefits  of  the  five 
multipliers  is  again  apparent  since  their  packings  are  consider¬ 
ably  better  across  dimensions  than  those  for  the  more  commonly 
used  multipliers. 

Knuth  (1981,  p.  102)  has  also  used  this  concept  of  packing  to 
rate  multipliers.  However,  his  approach  relates  to  packing 
spheres  in  the  dual  space  of  qQ/M . ^/M.  This  is  done 

A 

by  noting  that  in  addition  to  dfc(A,M)  being  the  maximal 
distance  between  neighboring  parallel  hyperplanes  in  the  space  of 

A 

Vfc,  the  quantity  l/Md^(A,M)  is  the  minimal  distance 

between  points  in  the  dual  space  of  q^/M,...,  q^  ^/M. 

* 

Therefore,  the  volume  of  a  sphere  with  radius  l/2dfc(A,M)  in 
the  dual  space  is 

ir  k/2 

W  (A.M)  = - - - 

r<k/2+l>  [2Mdk(A,M>] 

Now  observe  that  restrictions  (i)  and  (ii)  determine  that  the 
k  k  k-1 

hypercube  [1,1)  contains  exactly  2  M  k-dimens ional 

points  q/M.  In  particular,  the  exponent  k-1  instead  of  k  on  M  is 

due  to  restriction  (ii).  Therefore,  the  proportion  of  this  unit 
hypercube  packed  with  spheres  is 
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M.  =  2kHk'  V  (A,M) 
k  k 

k/2 

_ £ _ 

Rk/2  +  l)M{d*(A,M)}k 

which  is  the  measure  of  packing  in  the  dual  space.  This  quantity 
is  identical  with  the  figure  of  merit  suggested  by  Knuth  (1981, 
p.  101).  Note  that  because  of  the  lattice  structure  in  the  dual 
space  every  other  unit  hypercube  in  this  space  would  lead  to  the 
same  Wfc(A,M)  and  pfc(A,M) . 

Table  4  lists  p^(A,M)  for  the  multipliers  of  interest. 

Again  note  the  better  performance  of  the  top  five.  Knuth  remarks 
that  one  might  say  that  any  multiplier  for  which  p^(A,M)  >  .1 

k=2 . 6  passes  the  spectral  test  and  any  multiplier  for  which 

ufc(A,M)  >  1  k=2,...,6  passes  the  test  with  flying  colors.  By 
this  standard  the  top  five  multipliers  are  untouchable.  In  fact, 
since  S  (A.M)  >  .8  k=2,...,6  for  all  multipliers  in  the 

1  p  K. 

Appendix,  those  multipliers  have 


2.32  k=2 

3.03  k=3 

yk(A,M>  >  4.04  k=4 

4.88  k=5 

6.26  k-6 


indicating  that  all  meet  the  Knuth  criterion  and  dominate  all 
multipliers  listed  in  Table  1  of  Knuth  (1981,  pp.  102-103). 

Table  5  presents  bounds  on  discrepancy  computed  from  (25)  and 
(28)  and  reveals  several  interesting  results.  First,  note  that 
the  intervals  for  k=2  can  in  no  way  be  regarded  as  narrow. 


L 


Second,  the  top  five  multipliers  do  not  dominate  A=397204094  and 

630360016  unambiguously,  as  in  the  earlier  tables.  This  lack  of 

discrimination  on  the  part  of  the  lower  bounds  on  discrepancy  may 

be  due  to  the  fact  that  discrepancy  is  not  a  rotation  invariant 

measure.  That  is,  it  is  developed  along  the  lines  of  the 

classical  serial  test  in  Statistics  in  which  the  sides  of  the 

cells  are  parallel  to  the  coordinate  axes  and  hence  discrepancy 

detects  the  worst  case  with  regard  to  this  orientation  only.  By 
* 

contrast,  d^(A,M)  measures  the  worst  case  with  regard  to  all 

possible  orientations.  Although  one  can  argue  that  many  statis 

tical  testing  procedures  rely  exclusively  on  this  cartesian  pro 

duct  space  specification,  the  fact  that  our  study  reveals  so  many 

multipliers  that  perform  well  on  the  more  stringent  measure 
* 

dj_( A,M)  encourages  us  to  recommend  this  criterion  for 
general  use. 

As  mentioned  earlier  the  Appendix  contains  a  list  of  all 

multipliers  for  which  min  S  .  (A,M)  >  .80.  A  perusal 

2<k<  6  1 ’ 

of  this  list  reveals  six  multipliers  for  which  S  (A,M)  >  .80. 

3  *  It 

While  these  multipliers  do  not  rank  as  high  as  the  five  best  with 

regard  to  min  S  (A,M),  their  relatively  good  bivariate 
2<k<6  1,fc 

behavior  with  regard  to  ^(A,M>  and  S ^  ^_(A,M>  encourages  us 
to  examine  them  more  closely.  Table  6  shows  how  these  multi¬ 
pliers  perform  with  regard  to  lattice  packing  in  the  sample  space 


m 
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and  in  the  dual  space.  A  comparison  of  these  results  with  those 
in  Tables  2  and  3  makes  clear  that  these  multipliers  are  equally 
acceptable  with  regard  to  lattice  packing  considerations. 

Whether  or  not  some  other  justifiable  basis  exists  for  choosing 
these  multipliers  over  the  best  five  is  not  apparent  at  present. 

4 .  Empirical  Evaluations 

In  addition  to  evaluating  the  global  properties  of  a  multi¬ 
plier,  one  needs  to  consider  the  local  randomness  properties  of 
subsequences  of  moderate  length  that  a  generator  with  this  multi¬ 
plier  produces.  This  evaluation  is  usually  performed  by  statis¬ 
tically  testing  these  subsequences  to  detect  departures  from  ran¬ 
domness.  Fishman  and  Moore  (1982)  described  a  comprehensive 
battery  of  tests  for  this  purpose,  and  we  apply  the  same  battery 
here  to  test  the  five  best  multipliers. 

Recall  from  (2)  that  U  ,  U^,...  are  the  random  numbers 
normalized  to  (0,1).  Hypotheses  to  be  tested  include: 

H0:  {Uj ;  i=l,...,n)  is  a  sequence  of  i.i.d.  random 

variables . 

H]_:  {Uj;  i-l,...,n}  have  a  uniform  distribution  on 

(0,1) . 

h2:  (U2i-T>  u2i*  i=l,...,n/2  have  a  uniform 

distribution  on  the  unit  square. 

h3:  (u3i-2*  u3i-l>  u3i>  i*l. • * • . <n-2>/3  have  a 

uniform  distribution  on  the  unit  cube. 


H4:  Hq ,  Hj ,  H2  and  H3  hold  simultaneously. 
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For  each  multiplier  we  collected  100  consecutive  subsequences 

of  n=200,000  numbers.  For  each  subsequence  i  and  each  hypothesis 

3  a  test  statistic  was  computed.  Then  for  hypothesis  j, 

T. j , . . . »T^qq  j  were  subjected  to  the  battery  of  tests.  Let 

T„  have  continuous  cumulative  distribution  function  (c.d.f.) 

G.  under  hypothesis  j.  Then  G.(T..)  and  F .  .»1-G.(T. .)  are 
3  J  ij  i.J  3  ij 

distributed  uniformly  on  (0,1)  and  for  0<t<l 


F  .(t)  -  i  l  1  _(P.  .), 

n.3  n  (0,tl  ij 

where  I  denotes  the  indicator  function  on  the  set  B,  is  an 

O 

empirical  c.d.f.  If  H.  is  true 

3 

D  .  -  sup  |  F  .(t)-t  | 
n,j  t  n.j 

has  the  Kolmogonov -Smirnov  (K-S)  distribution, 


V  .  =  n  J  Irft  . .  (F  .  ( t) )dt 
n.j  Q  10, t]  n, j 

has  the  uniform  distribution  on  (0,1)  (Dwass  1958)  and  for  large  n 


A2  .  =  n  J1  UF„  .(t)-t]2/t(l-t)}dt 


n,j 


n,j 


has  a  distribution  given  by  Anderson  and  Darling  (1952,  1954)  and 

is  denoted  by  A-D.  The  quantity  D  .  measures  the  absolute 

n  >  3 

deviation  between  the  empirical  and  the  hypothesized  c.d.f.; 

V  .  measures  the  proportioon  of  F  .  that  lies  below  the 
n, j  r  n,j 

2 

hypothesized  c.d.f.;  and  A  .  is  a  weighted  measure  of  the 

ni3 

extent  of  deviation,  principally  in  the  tails,  of  the  empirical 
c.d.f. 
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Since  Fishman  and  Moore  (1982)  provided  complete  descriptions 

of  the  testing  of  here  we  merely  review  the  most 

essential  details.  In  particular  each  test  statistic  T.  .  was 

ij 

chosen  as  follows.  To  test  HQ  we  relied  on  a  comprehensive 

analysis  of  runs  up  and  runs-down  statistics.  For  we  chose 

12 

a  chi-squared  goodness-of-f it  statistic  with  2  =4096  cells. 

For  the  serial  test  statistic  was  used  for  nonoverlapping 

2- tuples  with  a  total  of  4096  cells  in  the  unit  square.  For 
Hj,  a  serial  test  statistic  was  used  for  nonoverlapping 

3- tuples  and  4096  cells  in  the  unit  cube. 

The  hypothesis  is  omnibus  in  character.  Recall  that 

P. .  =  1-G . (T. . )  i=l , . . . ,100  j=0,l . 3  and  set 

1 J  J  ^  J 

X. .  =  *_1(P. .) 

13  U 


where  ♦  1  is  the  inverse  of  the  unit  normal  distribution. 

Under  H.,  X..  has  the  unit  normal  distribtuion  and  X.., 

J  iJ  i0 

X^ . X_3  have  a  multinormal  distribution  function  »>.  Let 


X.  .  -  min  (X.  . . X.  _)  and  X.  =  max  (X.  „,..., X.  „). 

i,min  i,0  i,3  l .max  i,0  i,3 

Then  under  IT 


T.  .  =  l-*(-X.  .  ,  -  X.  .  ,  -  X.  .  ,  -  X.  .  ) 

1,4  i,min  i,min  l  ,mm  i,min 


and 


T.  A  =  1-HKX.  ,  X.  ,  X.  ,  X.  ) 

i,4  i,max  i,max  l .max  i ,max 

each  have  the  unit  r.brmal  distribution.  Since  T,  .  and  T.  ,  measure 

1,4  i,4 

how  likely  one  i6  to  encounter  values  as  extreme  as  X.  .  and 

l  ,mi  n 


X.  ,  they  provide  valuable  information  about  the  truth  of 

i  ,max 
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Ho . H3.  Accordingly  we  used  {T.  i=l . 100}  and 

*  »  ** 

{T.  i=l,...,100}  to  test  H,.  As  an  interim  result  a 

1,4  4 

test  of  the  multinormality  of  X^  q, . . . ,X.  3  was  also 
performed. 

Table  7  presents  the  P  values  for  and  the 

0  4 

multinormality  test  for  the  five  best  multipliers.  Although 
several  multipliers  show  some  small  P  values,  no  systematic 
rejection  occurs  across  the  K-S,  V  and  A-D  tests  and  across 
hypotheses.  If  one  feels  compelled  to  rank  the  multipliers,  one 
might  regard  A=950706376  as  first  and  A=1343714438  as  last. 
However,  we  emphasize  that  in  a  table  with  so  many  entries  some 
low  values  are  to  be  expected  when  all  hypotheses  are  true.  In 
summary  we  conclude  that,  in  addition  to  having  optimal  global 
properties,  the  five  multipliers  show  no  empirical  aberrations. 
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Table  1 

Linear  Congruential  Generators:  Zj  3  AZj  i  4  C(mod  M) 

Generated  Sequence 


Type 

M 

C 

A 

zo 

is  a  Permutation  of 

T 

1 

2b 

odd 

1  mod  4 

{0,1... • .Ml} 

{0.1 . M  -  1} 

2b 

2 

prime 

0 

primitive 
root  of  M 

U . M-l } 

{1 . M} 

M  1 

3a 

2B 

0 

5  (mod  8) 

1  (mod  4) 

{4j  +  l;  j=0,l,  .  . 

.  ,2b'2-  1} 

2B  2 

3b 

2b 

0 

5  (mod  8) 

3  (mod  4) 

{4jO;  j=0,l,.. 

.  ,  2B~  2-  1 } 

?B  2 

4a 

2b 

0 

3  ( mod  8 ) 

1  or  3  (mod  8) 

{8j4l  and  8343; 

3-0,1,.  . 

26-3-  1} 

2B  2 

4b 

2b 

0 

3  (mod  8) 

5  or  7  (mod  8) 

{ 8  j  -» 5  and  8  j  4  7  ; 

3  =0,1,  . 

,?B  3  1} 

2 

Source:  Jannson  (1966);  A,  C  and  Zq  guarantee  maximal  period  for  the  modulus  M=2B  with  B 
>  3. 


■ 
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Table  2 


Performance  Measures  for  Selected  Multipliers  in  Z[  i  AZj_x  (mod  M)a 

(M  =  23l-l) 


Dimension  (k) 


Multiplier 

A 

2 

3 

4 

5 

6 

742938285 

Si 

.8673 

.8607 

.8627 

.8320 

.8342 

S2 

.8362 

.6613 

.6618 

.6021 

.6075 

S3 

.8673 

.8751 

.8507 

.  7838 

.  7983 

950706376 

S1 

.8574 

.8985 

.8692 

.8337 

.8274 

s2 

.9211 

.8183 

.6555 

.6806 

.6822 

s3 

.8574 

.9093 

.8412 

.  7565 

.  7646 

1226874159 

S1 

.8411 

.8787 

.8255 

.8378 

.8441 

S2 

.8273 

.  7240 

.  7815 

.6492 

.6822 

S3 

.8411 

.8877 

.8468 

.7107 

.  7  743 

62089911 

S1 

.8930 

.8903 

.8575 

.8630 

.8249 

S2 

.  7169 

.  7537 

.7430 

.7153 

.6603 

s3 

.8930 

.8286 

.7712 

.8150 

.  7385 

1343714438 

S1 

.8237 

.8324 

.8245 

.8262 

.8255 

s2 

.8676 

.6404 

.6492 

.6702 

.7103 

S3 

.8237 

.7785 

.  7906 

.  7874 

.7747 

16807 

S1 

.3375 

.4412 

.  5752 

.  7361 

.6454 

s2 

.2565 

.3264 

.5714 

.6754 

.5888 

S3 

.3375 

.5404 

.6162 

.6187 

.5889 

397204094 

Si 

.5564 

.5748 

.6674 

.  7678 

.594  7 

S2 

.5966 

.5038 

.6239 

.6597 

.4206 

S3 

.5564 

.5543 

.  7302 

.  7849 

.6417 

630360016 

S1 

.8212 

.4317 

.7832 

.8021 

.5  700 

S2 

.8823 

.4373 

.6534 

.7173 

.504  7 

S3 

.8212 

.6354 

.6441 

.7983 

.5510 

7 

1000  Sx 

.1420 

4.882 

27.62 

78.13 

152.6 

1000  s2 

.1221 

3.413 

16.81 

41.19 

74.77 

1000  s3 

.1420 

.02650 

.02921 

.06746 

.2201 

\  s  V 

d*  (A,M) 

K1'*.  S2 

*  N*(A,M)/(k! 

M)1/k  and  S3 

*  c*(A,H> 

V  "l'k- 
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Table  3 

Packing  Measures  in  the  Sample  Space 
k/2  *  k 

<Dk  ( A,M)  =  ir*-  M[ck  (A,M)  ]  /T(k/2  +  1) 


Dimension  (k) 


Multiplier 

A 

2 

3 

4 

5 

6 

742938285 

2. 73 

3.97 

5.17 

4.40 

6.17 

950706376 

2.67 

4.45 

4.94 

3.69 

4.77 

1226874159 

2.57 

4 . 14 

5.07 

2.70 

5.14 

62089911 

2.89 

3.37 

5.17 

5.36 

3.87 

1343714438 

2.46 

2.80 

3.86 

4.51 

5.16 

16807 

.41 

.93 

.00 

1.35 

1.00 

397204094 

1.12 

1.01 

2.80 

4.44 

1.67 

630360016 

2.45 

1.52 

1.70 

4.83 

.67 

Upper 

Bound 

3.63 

5.92 

9.87 

14.89 

23.87 

-  39- 


Table  4 


Parking  Measures  in  the  Dual  Space 

k/2 


Plc(A,M)  = 


f(k/2-(l)  H[d*<A,M>]k 


(M=231-l) 


Dimension  (k) 


Multiplier 

A 

2 

3 

4 

5 

6 

742938285 

2.73 

3.78 

5.47 

5.94 

8.04 

950706376 

2.67 

4.30 

5.63 

6.00 

7.66 

1226874159 

2.57 

4.02 

4.58 

6.15 

8.63 

62089911 

2.14 

4.34 

4.23 

4.77 

7.99 

1343714438 

2.46 

3.42 

4.56 

5  .  73 

7.55 

16807 

.41 

.51 

1.08 

3.22 

1.73 

397204094 

1.12 

1.13 

1.96 

3.97 

1.06 

630360016 

2.45 

.48 

3.71 

4.94 

.82 

Upper 

Bound 

3.63 

5.92 

9.87 

14.89 

23.87 

ir 


-40- 


Table  5 


Bounds  on  Discrepancy 
Dimension  (k) 


Multiplier 

A 

2 

3 

4 

5 

742938285 

Lower® 

.1492 

.5970 

42.89 

42.89 

Upper*5 

3.446 

950706376 

Lower 

.2680 

1.072 

9.607 

10.08 

Upper 

3.725 

1226874159 

Lower 

1.967 

7.869 

7.869 

7.869 

Upper 

10.52 

62089911 

Lower 

.4236 

1.694 

1.694 

1.694 

Upper 

6.333 

1343714438 

Lower 

.2541 

1.016 

1.016 

1.016 

Upper 

3.772 

16807 

Lower 

1488 

5950 

5950 

5950 

Upper 

5952 

397204094 

Lower 

.4256 

1.702 

1.702 

1.702 

Upper 

4.517 

630360016 

Lower 

.1502 

.6008 

1.546 

1.546 

Upper 

2.980 

7 

Lower 

3571400 

14286000 

14286000 

14286000 

Upper  14286000 


aLower  Bound  =  108  x  1/min  (X_  t  |q,|). 

mU0  1 

h  «  P 

"Upper  Bound  =  10°  x  (2  +  £  a,)/T. 

U1 


6 


42.89 


10.08 


14.86 


4.328 


7.045 


5950 


28.61 


4.05  7 


14286000 
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Multiplier 

_ A 

809609776 

1567699476 

1294711786 

1554283637 

857010188 

1582405117 

Upper  Bound 


Table  6 


Packing  Measures 

for  Multipliers 

with 

Sl,k<A,M)> 

.  8  and 

S3(k(A.M)> 

.8  k= 

2 . 6 

Dimension 

(k) 

2 

3 

4 

5 

6 

Wfc(A,M) 

Wj.(  A,M) 

3.17 

3.17 

3.76 

4.23 

4.51 

4.55 

4.51 

5.07 

8.26 

6.71 

wk(A,M) 

Pk(A,M) 

2.88 

2.88 

3.66 

3.15 

4.98 

4.37 

6.55 

5.72 

9.96 

6.71 

wk(A,M) 

Pk(A,M) 

3.08 

3.08 

2.72 

4.73 

4.95 

4.73 

5.44 

4.17 

9.85 

5.65 

uik(A,M) 

Uk(A,M) 

2.56 

2.56 

3.71 

4.15 

4.71 

4.27 

6.08 

5.74 

7.79 

6.38 

wk(A,M) 

Uk(A,M) 

2.39 

2.39 

4.16 

4.20 

5.97 

6.39 

5.21 

5.95 

7.74 

5.08 

wk(A,M) 

Ufc(A.M) 

3.09 

3.09 

3.13 

4.25 

4.02 

5.24 

4.85 

4.88 

8.02 

5.78 

3.63 

5.92 

9.87 

14.89 

23.87 
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Table  7 


P  Values  for  Testing  Hypotheses 


Multi _ H4 


Multiplier 

A 

Test 

<1> 

H0 

(2) 

Hi 

(3) 

h2 

(4) 

h3 

(5) 

normality 

<6> 

min 

(7) 

max 

_ (8). 

742938285 

K-S 

.735 

.499 

.306 

.633 

.922 

.  776 

.802 

V 

.853 

.  012b 

.971 

.491 

.463 

.278 

.353 

A-D 

.408 

.231 

.406 

.  796 

.990 

.545 

.8/0 

950706376 

K-S 

.361 

.304 

.636 

.  766 

.163 

.244 

.529 

V 

.974 

.827 

.616 

.493 

.443 

.401 

.322 

A-D 

.269 

.254 

.497 

.629 

.173 

.279 

.417 

1226874159 

K-S 

.738 

.115 

.081® 

.903 

.151 

.220 

.532 

V 

.378 

.468 

.646 

.395 

.183 

.425 

.749 

A-D 

.442 

.083* 

.172 

.914 

.166 

.420 

.802 

62089911 

K-S 

.232 

.506 

.493 

.073a 

.578 

.121 

.132 

V 

.618 

.923 

.773 

.193 

.160 

.305 

.345 

A-D 

.328 

.457 

.539 

.139 

.377 

.151 

.  144 

1343714438 

K-S 

.771 

.068® 

.  024b 

.845 

.635 

.904 

.230 

V 

.849 

.440 

.158 

.781 

.577 

.365 

.404 

A-D 

.806 

.099® 

.041b 

.863 

.542 

.903 

.195 

a. 05  <  P  Value  <  0.1. 

b. 01  <  P  Value  <  .05. 
CP  Value  <  .01. 


M 
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Appendix 


A 

mn  Slfk 

min  S-,  . 
k  J,K 

A 

min  S,  , 
k  ]*k 

min  5^  . 
k  i,K 

7h2930285 

0 . 8  3  1 9 

0  .  7  >3  3  8 

1  76  00  2488  9 

0.0112 

0.  7  94  3 

5507063  7o 

0-82  7 4 

C.  75t>  5 

1  4422  73954 

0.81  11 

0. 71  10 

1228874  159 

u.  8  2!)  5 

0.  7  10  7 

9593874  18 

o. a  i  io 

0.77  4  5 

82dd9rJ  11 

O.  82*49 

0. 7385 

11 1  3  1  2  /  1 o  4 

0.01  08 

0.7726 

1343714438 

0.  82  3  8 

0.7747 

1448283090 

o. a  107 

0.7677 

2  04  95  1  3  9  12 

O'  .  8  2  3  2 

0.6  54  5 

231487336 

0.0107 

0.7020 

781259587 

U  -  d  2  1  2 

0. 709 9 

231487336 

0.81 07 

0.7820 

•♦0292c  3  80 

0. 8  2  04 

0. 7489 

4  0  3  6  3o  26  J 

0. d  102 

0. 794b 

18  108  <1  890 

U  .  8  1  9  8 

0.7852 

36  5. 3  70474 

0. 8098 

0.7375 

5  0  2  0  0  5  1 5  1 

O  .  198 

0.8930 

1  6  6  3  J  4  b  9  6  4 

0.809b 

0.71  13 

<♦  u  4  8  2  2  t  3  j 

0  .  o  1  9  1 

0.  7  3 8  8 

5o  4  t>  *  9  3  3 

0.  «09  5 

0.  7021 

1  4  o  0  3  o9Hr.H 

0-8  i  a u 

0.  7  34  5 

3  9  1  5  3  9  9  3  u 

0.8095 

0.7495 

.  3  2  4  4  "♦  0  4  1 4 

i;  •  h  1  o  4 

0. 7271 

621 -89603 

0.8093 

0. 7b7o 

1 9302 51 322 

u.ll  e  2 

0.7159 

18978  3 o 7 2 2 

0.  8092 

0.7b  Io 

800  218  25  1 

0  .  f>  1  o  2 

0 • 7  3b  6 

209  72044  3 

0.80  92 

0.7582 

1  5  /  59t.  5  843 

0.8181 

o.  7242 

1n51 1 32469 

0.8090 

.7.7605 

1 1  0  04  >4401 

o . a  1  70 

0.0828 

10  3t  4(397  97 

0 . 509u 

0.7361 

1  (>4  72  74  979 

1/  m  O  1 

0.  7  124 

1 0940  0..H  55 

0. 80 8c 

0.7044 

o  2  2  5  2  5  8  d 

0.8188 

0.  7  35  4 

9583732  00 

0.0088 

0.7173 

1  9  0  4  5  o  5  5  2  9 

d  -  6 1  6  O 

0.7577 

1 8  8  2  4  u  *  2  1  d 

0. 808  / 

0. 7956 

10  12  14  3948 

u  .  ft  1 6  4 

0. 7470 

19019 18329 

0.0  08  7 

0.7586 

1  7  j  4  0 8  04  0 0 

6  .  f:  155 

0. 7455 

1482800924 

0 . d  0d4 

0.7763 

1  5  o  C  c  5  0  o  3  o 

6  .  ft  1  J‘4 

0.7  72  J 

1 6  2  9  2  d  o  0  5  1 

0  .  d  0  7  >' 

0.  74  30 

1022204 322 

6  .  4  1  5*4 

0.7076 

18734^6901 

0. 8u 78 

0.6724 

30  0  1 0  ri  o  1 

0 . 0  1  j  2 

0.  744  1 

1J 9 46 33640 

0 . 8075 

0.7039 

1  1  {  7  d  4 1 4  5  3 

0 . 8  1  -1  0 

0.7312 

1  o  9  1  9  1  j  .01 

0  .  u  0  7  5 

0.7119 

5  3  17  4 9  2 4.5 

J  .  8  1  <4  8 

0.  7  17  9 

155.  7 9  ■  c 2 

0.8075 

0.6  7  H 

1  40  25  5  1 8  14 

0.0147 

0.727 7 

1  4095  5jo<.>  7 

0 .  (J  0  7  3 

0.7992 

5  c  8  7  4  5  7  5  7 

0  .  o  1  *4  J 

0. 7587 

2 1 1 1 J On  7  2  1 

u .  30  7  3 

0. 7198 

1  o  i>  7  4  0  j  9  16 

O  .  '■ '  1  4  4 

0. 7545 

1  J  37  2  39  1  3  9 

0.8072 

0.  7897 

14344725 91 

0 . 8  1  -42 

0.751 7 

1257701541 

0.8072 

0.7358 

1  5428  73971 

9  .  !3 1  *4  2 

0.7938 

1  06  1 0 2  3  7 9ft 

0.8072 

0.  7087 

o  2 1 506  5 30 

u  .  o  1  *♦  1 

0. 7  158 

65994/2/0 

0.807  1 

0. 668  9 

473911478 

0  .  c 1  3  9 

0. 7548 

1  4  72-J  0/  7 1  -  o 

0.8071 

0. 7432 

2  1103  r.2  50t> 

6.  31  3  9 

0.7  78.3 

1 7  09  3  5 *♦  -*  o  2 

0  .  33  0  o  9 

0.7457 

1  5  0  O  O  3  8  4  O 

0  .  o 1  3d 

0.7012 

1  43  73  552  12 

0. 8069 

0.7240 

1  J  In  4844(3 

0  .  o  1  3  8 

0.  7  74  0 

21121 59807 

0.80b 9 

0.7122 

1  1  1  4950053 

8.(31  3  3 

0.75OH 

1  6  1  0  3  >0  8  1  0 

0  .  0  0  6  8 

0.70  29 

17 1, 80  533  39 4 

0  .  0  4  3  0 

0.  7  50  9 

1  36232  564*. 

0. 0068 

0.6  36  9 

51  3482507 

0  .  ft  1 2  7 

0. 78uJ 

1  52  8 1  oo  t.  52 

6 . 0  06  d 

0 . 7  7  /  8 

5  134c2  or- 7 

0.8  127 

0.7803 

o  4  4  5  1  3  4  7 

0.  8  06  7 

0. 785b 

1o2o240()4  5 

0.0127 

0.  7  30  8 

1 640  0  113  1a 

0. 8063 

6.7232 

2  0  99  4  8  5  7  54 

j .  o  1 2  7 

0.  740  8 

1  2  6  7  2  0 1  1  7  0 

0. 6062 

0.7771 

K  o  2 1  1  J  d  1  8 

J  .  8  J  2  / 

0.8294 

009009776 

0.0  0b  1 

0.0222 

J  3  4 \t  3  1  1  48 

u. ft  125 

0. 6  d4  9 

29. 357876 

6. 006  1 

0. 7322 

*40  42  c8  709 

0.-3  124 

0 .  7  20 1. 

1022131410 

0.8 06 1 

0.7509 

.  5  7  2  o  C  J  3  9 

0.0  1  2*4 

0. 7 Jon 

18366 24202 

0.8061 

0.7595 

100o097‘4OJ 

o  .  (3  1  2  1 

0.7780 

67  2  :5  3o  7  1  7 

0. 8060 

0.  7532 

1  5  9  j  4  9  2  7  5  / 

0.8121 

0. 7484 

1  2  92  U  6  8  b  4  8 

0. 8059 

0. 66  /  3 

48  4  0  2  c  2  80 

O.i.  11  5 

0.7029 

90  5  146  4 0  4 

0  .  d  0  5  9 

0. 7546 

•f* 

The  remaining 

207  multipliers  can 

be  computed  as 

fol lows : 

Set  B  = 

for  each  multiplier  A  find  the  smallest  integer  1  such  that  A  1$' (mod  M) . 
Then  the  multiplier  A*  E  B^-l"I(mod  M)  has  the  same  properties  as  A. 


Appendix  (continued) 


A 

min  S,  , 
k  1>k 

min  S.  . 
k  3’k 

149  38 J4fc01 

.  t  i  0  b  7 

0.  6  70  b 

1  0  3  75ub  900 

•; .  oo bd 

0. 7407 

7  4  37  2 4  n h 

j. .0  bo 

0.  7o j7 

1  bo  90  899  i  1 

b  .  o  0  9  7 

0. 72o  4 

1  5 1 j  7 «  5  S  4  7  6 

o . j 7 

0  .  8  4  2  o 

194730b9.<7 

o  .  c  0  b  3 

0 .  7  1  u  4 

10  / t  5  J  l  C  9  / 

J  .  •  ■  J  b  2 

0. 7b0  J 

1 9b m 1 1 7J  / 

0  .  ( •  0  \>  . 

G. b  03  4 

b  2  8  4 '» / 1 4 f ) 

J  .  oo  b  1 

0.  7  54  0 

1  0  4  'J  () 4 5  3  9  3 

0 . 0047 

0.  727.: 

7  >3  (,  8  2  4  *  .3  _> 

U  .  ■  j  .J  J 

0. 770' 

5  5  b  5  >  0  h  2  4 

\j  .  O  0  t  7 

0 . 7  3*0 

37b. 1 2 jO 

0  .  c.  )4  / 

0 • 74  o  2 

1 4b 7b 13431 

u  .  j  ./4  7 

0. /Vo  6 

J  d  5  7  d  7  4  b 9 

0  .  O  0  4  0 

0.  7b7o 

15t>73l65l2 

0  .  ‘  0  4  Q 

0  .  7  bo  o 

l)  3  C  9  j  5  3  4 1 

0 .  o04  4 

0.7  77  0 

1  56  88  13  46b 

0  .  u  0  4  4 

0. 7850 

103  55  152  17 

0 . 3  0  4  3 

0.  7bV0 

3  9  4  4  J  4  b 

0.0043 

0.6732 

189  133047  3 

0  .  hj  4  3 

0.  705B 

169/4  12292 

o. 3043 

0. 7112 

7  54ocCJ39 

0  .  1 0  4  3 

0.744  7 

1 9  7 12048  12 

0  .  o  0  4  3 

0.7733 

166  88  >  7  7  9  o 

o  .  tfO  42 

0. bbj  d 

1  J  7  1  o  4  1 t  J  4 

0.8  )4  0 

0.744b 

1  1  1  74  .35  554 

0.  oO  4  0 

o.  7  24  J 

bob  1  70  t  ti2 

0  .  c  0  4  0 

0. 7292 

9  2  7 4  072  j  4 

0 . 0  0  4  0 

7.7147 

1  4()uob00o  7 

0  .  (>>}  3  7 

0.1220 

2  3  5  7  It  S  7  7 

0 . 00 J  7 

0.731  J 

1  4  9  2  8  9  fc  2  b 

0 . t03 0 

0.7  02  o 

166v>3  7 1;  127 

o .  in' *2 

0.  7ri5  1 

15172061)7 

O.oO  3d 

0 .  r>  32  7 

12obdH1  C  1 .! 

o  .  d  0  j  7 

)  .  7  14  6 

7  0  7 1  bo  2  77 

■o  ,t.O  <  7 

0.7017 

1  H  ti  70  7  5  7  34 

o . oo 3  7 

0.6714 

b7b JO04O4 

0  .  d  0  j  7 

J .  7  1 6  2 

b  37 1 4t  2i>B 

o  .  d  0  J  7 

0. 7  50  3 

It. 04  1c7  177 

o  .  oo.3 n 

0. 701 3 

JOd^l 20220 

o. d03 b 

0 . 7624 

.1705447  24 

U.  d  0  J  b 

0. 73 7b 

2  0  4  4  7  2  4  5  7  1 

0. d03 o 

0. 6  9)38 

71 >1 0  C  7  u  7 

J.odo b 

0.6  07  9 

1o  3  74  14  1*» 

u.  oO  jb 

0  .  7  d  6  6 

13  3c.  1.2  4  10 

O.otl  3  3 

0.  7246 
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